N79-30048 

18 

SUPERCRITICAL TESTS OF A SELF-OPTIMIZING, 
VARIABLE-CAMBER WIND TUNNEL MODEL* 

Ely S. Levinsky 

General Dynamics Convair Division 

Richard L. Palko 
Sverdrup/ ARO, Inc. 


SUMMARY 

A testing procedure has been used in the 16-foot Transonic Propulsion Wind Tunnel (PWT-16T) 
at the Arnold Engineering Development Center (AEDC) which leads to optimum wing airfoil 
sections without stopping the tunnel for model changes. Being experimental, the optimum 
shapes obtained incorporate various three-dimensional and nonlinear viscous and transonic 
effects presently not included in anal 3 rtical optimization methods. The present method is a 
closed-loop, computer-controlled, interactive procedure and employs a Self-Optimizing Flexible 
Technology (SOFT) wing semispan model that conformally adapts the airfoil section at two 
spanwise control stations to maximize or minimize various prescribed merit functions (e.g., 
miniimim drag) subject to both equaliti^ and inequality constraints (e.g., fixed lift, maximum 
spanwise differential deflection). The model, which employed twelve independent hydraulic ac- 
tuator systems (nine functioned) and flexible skins, was also used for conventional testiiut (Ref. 
1). Although six of seven optimizations attempted were at least partially convergent (several of 
which are shown herein), further improvements in model skin smoothness and hydraulic 
reliability are required to make the technique fu^ operational. 


INTRODUCTION 

Although considerable interest has been generated at this Advanced Technology Airfoil 
Research conference and in the past in the anatytical design of optimum airfoils (e.g.. Ref. 2-4), 
such methods are somewhat approximate, because of the inability of current-generation com- 
puters to evaluate the flow field with sufficient accuracy and rapidity. This is especially true 
under conditions of high lift or supercritical transonic flow, where boundary layer transition, 
separation and possible reattachment, and shock wave interaction are among the highly non- 
linear flow phenomena still beyond analytical solution. In addition, the large number of flow 
field solutions required for each optimization may preclude analytical optimization of three- 
dimensional wings, aven with the most advanced digital computers. 

The present SOFT wing optimization procedure differs from the analytical methods in that 
the wind tunnel is used as an analog computer, in corjunction with a flexible and controllable 
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model, to rapidly generate aerod 3 mamic data, instead of using a digital computer for this pur- 
pose, In both cases, the data are then supplied to a computer, which controls and optimizes wing 
shape by application of nonlinear programming techniques. The SOFT wing procedure will, 
however, be limited by the extent of wing articulation and degree of shape control that may be 
‘obtained with a wind ttmnel model. 

Under a previous contract for the U.S. Navy (ONR), General Dynamics Convair Division 
fabricated and tested a two-dimensional flexible airfoil model that could be optimized in the 
tunnel in a manner similar to that used for the present SOFT wing model. The 2-D model was 
tested under both low-speed and transonic conditions (Ref. 5), and employed five pairs of hy- 
draulic actuators to vary leading edge radius and camber. Airfoil shape and angla of attack were 
con^lled “on line” to minimize drag, subject to constraints, using the gradient projection op- 
timization algorithm. Satisfactory convergence was found, even with flow separation; however, 
becatm of the intermittent blow-down type operation of the transonic test facility, a sequence of 
individual runs was required to complete a single optimization problem. 

The present 3-D SOFT wing model was tested in PWT-16T at AEDC during the summer of 
1977. Test Mach numbers ranged from 0.6 to 0.925. The SOFT wing model employed twelve in- 
dependent hydraulic actuator systems similar in design to those in the 2-D model. The com- 
puter-controlled testing and optimization techniques were also similar, but the number of con- 
trol channels was increased to handle the larger number of independent control variables. 
Because of the continuous operation of Tunnel 16T, continuous computer control of the wing 
could be maintained during an optimization, and under ideal circumstances, a single run in the 
wind tunnel was sufficient to arrive at an optimum wing shape for a prescribed test con ditio n , 

The present paper briefly describee the SOFT wing model, closed-loop testing technique 
and optimization procedure, summarizes several of the optimization problems attempted, and 
reviews the causes of various difficulties. 


SYMBOLS 

Ay Actuator positions in counts (0 < Ay ^ 1000) 
Cq Drag coefficient 

Lift coefficient 

F Objective function being minimized 

g Vector of constraints 

Kj Stepsize of i^^ simultaneous mode point 

M Mach number 

N Iteration number 

R Restoration function 

S,Sr Vector search directions 

a Angle of attack 

A Lagrange multiplier vector 

<l> Merit function to be minimized 

9 Independent variables (transformed) 

V Gradient vector 






-/ .. ^ 




Superscript 

T 


Transposed matrix 

Best K; also uncorrected for weight tares 


WIND TUNNEL MODEL 

existing horizontal tea flow-through naceUe, an 

edge sweep angle 26*. The wing actuators sho»m u * leading 

perform the following variations: ’ schematically in Figure 1, are designed to 


^11 

Ai2 

Ai3 

Ai4 

Ai6 
Ai6 
A21-A 


■26 


Nose radius, inboard station 

N«e deflection about 16* chord line, inboard «ation 
No» deflection about 26* chord line, inboard station 
Upper surface humping, inboard station 

Trailtog edge deflection about 66* chord line, inboard station 
Traihn* ^ deflection about 80* chord Hne, inboard station 
ftixiB as 11 ‘Ajg, 6xc6pt at outboard stations 


All systems except for the Ae . 4 sa«j a a . 

diving a total often degrees if frSom todudS‘X'1^^^ «>' *“1. 

ported from the win?tUJ!^l sWew^l*%e mo^T^h^^®/ reflection plane sup- 
tested in Figure 3. assembled and installed as 

CLOSED-LOOP TEST PROCEDURE 

were'tostaZl^i!i!fid‘^^n/^/".gi ^ the time that the skins 

ment. of the control system and »m”^ B^^^l^ll ‘k* •»«*' 

through the PWT DEC System 10 comnuter wi. f ‘ Position changes were made 

venrional runs were made either by steppiig the artuTtow “k ** Con- 

positions at fixed a and M or bv hoi5?«!?*k through a prescribed sequence of 

prescribed schedule. ' ^ ®nd stepping « thZ^ a 

tuns «rSctnS^trr«SSr^^^ ^«ween 

ate<'tl^*-Srr^ in that it goner- 
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OPTIMIZATION PROGRAM 

A flow chart of the optimization procedure, which was based on the gradient projection op- 
timization algorithm, is presented in Figure 6. The optimization program operates in two dis- 
tinct modes: “incremental” and “simultaneous.” Only the incremental mode is used during the 
initial iteration (and restarts), during which each active actuator and angle of attack is per- 
turbed individually to generate gradient vectors V.0 and V g of the merit function and active 
constraints, respectively. The number of incremental mode points per iteration depends on the 
number of active actuators and on the number of times the perturbations are repeated to im- 
prove accuracy (termed “cycling”). With 9 active actuators, there were 11 incremental mode 
points per cycle. 

After the incremental mode, the vector directions Sg or S, for either restoring the con- 
straints g » 0 or minimizing the objective function F during the next iteration, depending upon 
whether any constraints were violated at the nominal incremental mode point, are calculated by 
the gradient projection algorithm, which gives in matrix form: 

F--0+gA 

S — VF 

Sr- v«(Vg’’v«)“‘g 

and 

A--(vg'*’v^*‘ 7g^V0 

The next iteration begins with the simultaneous mode, during which all a«.tive actuators and a 
are advanced together in the direction S or Sr, as obtained from the previous iteration, through 
a sequence of up to 11 test points of step size Kj. The sequence is aborted if any of the con- 
straints are violated by more than the prescribed tolerance. Upon completion of the stepping 
the computer selects the “best” of the simulteneous mode points K* (either minimum R - gg^ 
or minimum F), aiid then resets the wing shape to that configuration. The incremental mode is 
then repeated in prepairation for the next iteration. 

CONVENTIONAL RUNS 

Only the conventional run data comparing the SOFT wing T -1 shape with a comparable Mlid- 
wing model equivalent to the V/52 theoretical shape shown in Figure 2 will be presented. Figure 
6 shows that wing shape T-1 experienced significantly higher drags, particularly at the higher 
Mach numbers, even though T-1 was set to match the W52 theoretical shape as closely as possi- 
ble. The higher drag was most likely caused by the leading edge sliding ekin joint (lap) of ap- 
proximately 0.13 cm thickness that protruded out of contour on the upper surface, by a bulge 
that protruded out of contour in the outboard region of the aft lower surface, and by deviations 
in nose shape. The deviations in Figure 2 between the pre- and post-test T -1 shapes may account 
for the scatter band in Figure 6, since the flow at transonic conditions is highly sensitive to small 
changes in airfoil shape (Ref. 2). It is planned to improve the SOFT wing model for subsequent 
tests to approach solid model drag levels. 

OPTIMIZATION RUNS 

The optimization problems are summarized in Table 2 , which lists Mach number, merit func- 
tion, active constraints, convergence assessment, iteration number N considered optimum, and 
the percentage reduction in merit function. Each optimization problem was initiated with the 
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wing set at the T -1 target shape; the percentage of decrease is with respect to the T -1 shape 
tested dunng the same (or closest) run sequence. 

Weight tare corrections were inadvertently omitted during the on-line data reduction used 
to st^r the direction of mng optimization Consequently, evaluations of convergence were 
based upon uncorr^ted coefficients and C^, rather than on Cl and Co- This error biased the 
optimwation procedure in favor of increasing a at the expense of aft camber 

Only optimization 4-2 (completed as' Problem 30) will be discussed in detail. It 

tn~ arv? Trj ^ Subject to various inequality constraints on the actua- 
tore at Mach 0 86 . As shown m Figure 7. Problem 4-2 was continued for six iterations, after 
which a premature shutdown was experienced, due to lubrication difficulties with a main com- 
pressor beanng. The optimization was then continued as Problem 30, without reinitialization. 

J/ .r ® *^P*^*®®* representation of the convergence process. All constraints are 
Mtisfied for the iterations represented by solid symbols, while the open symbols signify that at 
least one constraint has been violated. Iteration 17 is readily seen to have the lowest value of Cn 
satisfied; It was, therefore, selected as optimum. The gap between Iterations 

fi ^ ** ^ ^ shutting down during Iteration 

6 . The double values for Cp during Iterations 12 and 14 and the poor convergence displayed dur- 
ing this ponod were the result of excessive shifts in the Ajs actuator position 

dependent data functions during 
Iteration N « 11 from Problem 30 are shown in Figure 8 . The iteration, during which Cn is 
minimized, contains 11 simultaneous and 22 incremental mode points cycled twice. The precise 

P^«*twns and the good repeatability of the aerodynamic data during the two 
incremental mode cycles are apparent. 

f Co plane is shown in Figure 9. Also shown are the drag polars for 

thS optimization run, and that made with the op- 

timum wing 0-30 (Iteration 17). The advantage of the 0-30 wing shape over T -1 at the desim 
Cl is clearly apparent. Unfortunately, this same reduction was not found in Cq because of the 
tare to previously. A correction for tare effects (made after the test for iteration 

“ ^ showed that increased trailing edge deflections, decreased angles of attack, and 

decrMsed values of Cp would have been obtained with the tare correction included (Figure 10 ) 

A q 18 for problems 3-2, 41-2, and 

4-3 of Table 2 . Problem 3-2 is similar to 4-2 and 30. except for reduced articulation (only the 

Ai 2 » a, 5 , A 22 , and A 25 actuators were varied, subject to A ,2 » Aoo and A,r - AoO Fieure 11 

itrThlTeducorH^^^^ ‘*'® (N - 20) was now much 

-IV. u ^ ^ articulation is apparent in Figure 12 (iteration N - 18), which also shows that 

the simultaneous mode was aborted after the ninth point of this iteration because the constraint 

o ; o ‘**® reduction for the optimum wing shape 

designated 0-3.2 persisted over a wide range of Cj,. ^ 

*^oWem41-2dealtwithminimizingC5atC£;'- 0.26withfullmodelarticulation.Excellent 

Iteration N - 29 was selected as optimum (Figure 14). Iteration 
N - 21 (Fi^re 15) is an example of a restoration. Note that all active constraints Ct and lA - 

‘**® “Ode- The convergence if Problem 

41-2 in the Cl, Cp plane is shown in Figure 16. 

Optimization Problem 4-3 was similar to 4-2 and 30, except that M - 0.90 instead of 0.85 

ThUTi®"'® *j-^**“ “®® N - 10 selected as optimum. 

Jeen^i^FIgJ*"^?* ******* ^ ***°'^“ “ considerable improvement over the T -1 wing. 
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The resultant optimum wings are compared with T-1 in Figure 19. Although the weight 
tare biased the optimum shapes toward decreased aft caml^r, the increased leading edge droop 
of the optimum ahflp eg, especially outboard, may well signify a shape modification that leads to 
lower drag at moderate to high Cl* On the other hand, Problems 42 and 10 (Table 2) showed 
m argiwQi ami no convergence, respectively. Both problems were attempted at high Cl and 
showed evidence of reduced control precision, possibly due to the high airloads on the actuators 
and/or buffeting. Further testing with an improved SOFT wing model is planned to oveifome 
some of these difficulties and make experimental wing optimization an operational tool for the 
aircraft designer. 

CONCLUSIONS 

A first attempt has been raade at developing an experimental 3-D wing optimization pnx^ure 
involving a computer-controlled SOFT wing wind tunnel model. Although six of seven optmaiza- 
tion problems attempted were at least partially convergent, difficulties were uncovered with 
deviations in airfoil contour, inadequate control precision at high Cr^, and hydraulic system 
reliability, which increased <h»g levels and slowed or prevented convergence. Nevertheless, the 
SOFT wing technique appears to promise a means of generating optimum airfoil and wing 
fthw pftg, which include all aerodynamic nonlinearities and viscous effects, with a saving in test 
fimo ami which might not be found by conventional testing or by numerical optimization. Of 
course, any optimum wings will be limited by model articulation, and will reflect effects of tim- 
nel wall interference, test Reynolds number and model deformation under airload. Plans call for 
improving the model and performing further tests to make the technique more operational. 
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TABLE 2.- SUMMARY OF OPTIMIZATION 





PROBLEMS: 

AEDC 457 



Problem Mach 

Merit 


Convergence 

Beat 

% C{| Decreaae 

No. 

No. 

Function 

Active Conitraint^ 

Propertiee 

N 

from SUndard Wing 

30 

0.66 

MlnCo 

C|, - 0.5, A Twiet ^ 4" 

Good 

17 

18% 

3-2 

0.85 

MinCQ* 

Cl = 0.5, “Reduced” 
Articulation, A Twiet ' 0° 

Excellent 

20 

8% 

41-2 

0.86 

MlnCo* 

Cl^ 0.25,ATwirt» 0 ± 2« 

Excellent 

29 

17% 

4-3 

0.90 

MinCo* 

Cl ' 0.50, A Twiet ^ 3** 

Excellent 

10 

14% 

42 

0.85 

MinCp* 

Cl 0.70 A Twiet == 0 -t 2* 

Very marginal 
except for laet 
4 iteratlone 

33 

9% 

6-2 

0.85 

MinCp* 

^ 0.50 

A Twiet ^ 3® 

Excellent up to 
N- 8 

Poor for N > 9 

6 

18% 

10 

0.90 

MaxC|; 

Cp** 0.ll.ATwiet^ 3* 

Poor 

None 

X 


Cd* " Co-0.04Mn(d«f) 
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Figure 2.- T-1 actual and theoretical wing shape. 
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Flow chart of mathematical optimization procedure. 



6.- SOFT wing and solid wing drag variations. 
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7.- Optimization summary. Problems 4-2 and 30. 
Cjj at Cj^ - 0.5; Mach 0.85. 
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Figure 9.- Problems 4-2 and 30 in C*, C* plane. 
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Figure 10.- Tare corrections for iteration 4 of problem 4-2. 
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Figure 11.- Optimization summary. Problem 3-2. Minimum C at C =0.5; 

limited articulation; Mach 0.85. ® ^ * * 
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Figure 12.- Minimization of C_. Iteration N - 18 of problem 3-2. 
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Figure 14.- Optimization summary. Problem 41-2. Minimum C at C 

Mach 0.85. D 
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Flgure 15.- Restoration if equality constraints on C. . Iteration 

N « 21 of problem 41-2, 



Figure 16.- Problem 41-2 In C*, C* plane. 
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Figure 18.- Problem 4-3 in C*, C* plane. 
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Figure 19.- Optimum wing shapes 



